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Advances in neutrino transport allow realistic neutrino interactions to be incorporated into su- 
pernova simulations. We add tensor couplings to relativistic RPA calculations of neutrino opacities. 
Our results reproduce free-space neutrino-nucleon cross sections at low density, including weak 
magnetism and recoil corrections. In addition, our opacities are thermodynamically consistent with 
relativistic mean field equations of state. We find antineutrino mean free paths that are consider- 
ably larger then those for neutrinos. This difference depends little on density. In a supernova, this 
difference could lead to an average energy of £> M that is larger than that for by an amount that 
is comparable to the energy difference between and 

I. INTRODUCTION 

Core collapse supernovae are dominated by neutrinos, therefore many supernova properties may depend on the 
neutrino-nucleon interactions. Recent advances in neutrino transport allow neutrino interactions to be accurately 
incorporated in simulations [1-3] . It is important to systematically improve the treatment of neutrino-nucleon inter- 
actions in order to take full advantage of the high quality neutrino transport and insure that modern simulations are 
as meaningful as possible. 

There are fundamental differences between neutrino- and antineutrino-nucleon interactions related to parity and 
charge conjugation violation. These differences were incorporated in refs. [4,5] by including weak magnetism correc- 
tions. These corrections are important for both charged and neutral currents. Weak magnetism has been included in 
supernova simulations with an average correction factor [6,7] and leads to a lower average energy for compared to 
Vfi [8]. The average neutrino energy for P M , v T was found to be over one MeV larger then that for v^, v T [8]. This 
significant difference is comparable to the approximately 1.5 MeV difference in average energy of v x (all four P„, D T1 v„ 
and v r combined) compared to v e found in a modern simulation that includes a variety of neutral current interactions 
such as nucleon bremsstrahlung and v e , v e annihilation [9]. 

Many-body corrections to neutrino interactions may also be important in the dense protoneutron star. Nonrelativis- 
tic calculations [10] find that correlations reduce the neutrino opacity at high densities, while relativistic calculations 
[11] find a similar but somewhat smaller reduction. Unfortunately these corrections depend on the model employed for 
the strong interactions. Furthermore, the model of the strong interactions used for the equation of state (EOS) needs 
to be consistent with the neutrino opacities to insure thermodynamic consistency. For example, uniform neutron rich 
matter is unstable at low densities to a phase transition to a nonuniform crust phase. This should show up as an 
enhancement of neutrino opacities as neutrinos excite low-energy density oscillation modes. 

We believe the minimum requirements for accurate neutrino opacities are the correct low density limit and con- 
sistency with the EOS. First, the neutrino interactions, at low densities, should agree with the exact free-space 
neutrino-nucleon cross sections including weak magnetism and recoil corrections. This will insure that fundamental 
differences between neutrino and antineutrino interactions are included. Second the opacities should be consistent 
with the EOS employed in the simulation. This will insure that the opacities correctly reflect possible phase transitions 
in the EOS. To our knowledge, these very minimum requirements have never been satisfied in a supernova simulation. 

Most supernova simulations use the Lattimer Swesty EOS [12] based on a nonrelativistic model for the strong 
interactions. It is important to study the sensitivity of simulations to changes in the EOS. In addition, we wish 
to include weak magnetism using a relativistic formalism. Therefore we consider the EOS [13] based on the TM1 
relativistic mean field interaction [14,15]. We calculate consistent neutrino opacities for this EOS in a relativistic 
RPA approximation. Our starting point is the relativistic RPA calculations of Horowitz and Wehrberger at zero [17] 
and finite temperature [18]. Previous relativistic RPA calculations with the TM1 interaction did not include weak 
magnetism [16]. Some additional relativistic calculations of neutrino opacities include [19-22]. 
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Section II presents our formalism for the neutrino cross section including weak magnetism corrections. Results for 
mu and tau neutrino transport mean free paths are presented in section III as a function of density and temperature. 
We also discuss differences between neutrinos and antineutrinos. We focus on the low density region between the 
neutrino sphere and normal nuclear density. These results are discussed and we conclude in Section IV. 



II. FORMALISM 



In this section we calculate neutrino scattering from neutron rich matter. We assume only neutral current interac- 
tions which is appropriate for /j, and r neutrinos and antineutrinos. We use a relativistic field-theoretical model for 
describing a system composed of protons, neutrons and electrons. 

In this work the TM1 model [14] is used where the nucleons interact via the exchange of a,ui and p mesons. 
Non-linear scalar and vector self-interactions are included. The TM1 nucleon- meson couplings are known to provide 
excellent descriptions of the ground states of heavy nuclei including unstable nuclei [14]. The interaction Lagrangian, 
in our conventions is, 



Ant = "0 



V„ 



r 



K 



a) 



The model contains an isodoublet nucleon field {jp) interacting via the exchange of two isoscalar mesons, the scalar 
sigma (<j)) and the vector omega (V^), one isovector meson, the rho (b M ), and the photon (A^). The parameters are 
listed in Table I. 

The differential cross section per nucleon for neutral current ^-nucleon scattering is [17], 
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where is the baryonic density and E (E') are the incoming (outgoing) neutrino energies. The energy transfer 
is qo = E — E' . In terms of the incoming neutrino four momentum fc M = (E, k) and the four momentum transfer 
<Zm = (qo,<l), the lepton tensor L a p is given by, 

L af3 = %[2k a k + (k ■ q)g afl - {k a q f3 + q a k fi ) T ie a ^ u k^q u ] (3) 

with a minus sign for neutrinos and a plus sign for antineutrinos. 



A. Hartree Response 



The retarded polarization tensor 11^, is a function of the particles chemical potential pi, temperature T, and the 
kincmatical variables qo, and |<f|. For a plasma of protons i = p, neutrons i = n, and electrons i = e it is, 

Imn^=tanh(|,) £ Im U% (4) 

where IT^ is the time ordered or causal polarization and is given by, 

H$ = -< / [T(G t (p + q)r a G t (p)r )]. (5) 

The Green's functions Gi(p) describe the propagation of baryons (i = p or n) or electrons (i = e) at finite density 
and temperature in the several approximations considered. 



TABLE I. Parameters for the TM1 interaction using the conventions of Eq. (1) 
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In the Hartree or mean field approximation d is [17,18], 

G,(p) = (fi + M*){ } + ^rS( P * E* p )n;}, (6) 

Pi lvl i lt pi 

with E* pl = (p 2 + M* 2 ) 1 / 2 , p* Qi ^pa-V, and p* = p. The mean fields are, 

V n = g v V - g P b /2, (7) 
V P = g v V + g p b Q /2, (8) 

V e = 0, (9) 

with Vq the Omega (vector) and b the Rho mean fields and g v , g p the respective coupling constants. The effective 
mass is M* = M* = M — g s <p for i — p,n and M* = M e . The scalar mean field is 4>q and the scalar coupling is g s . 
Finally n p ^ gives the occupation of a state with momentum p, 



{l + exp[(£& + V--/i i )/T]} . (10) 
The weak interaction vertex is F % . It includes vector, axial, and tensor (or weak magnetism) contributions, 

rj,(9) = FKQ 2 )^ + iFi{Q*)o^ - G\(Q 2 )^ 5 . (Q 2 = q 2 g 2 ) (11) 

The form factors i^, Fj, are given in Table II. For astrophysical purposes we can neglect the Q 2 dependence since 
the medium is being probed at Q 2 « 0. 

In this paper we will be interested in the Hartree (H) or Mean Field approximation and in the Random Phase 
Approximation (RPA) where correlations are included. The Hartree and RPA response functions which characterize 
the hot and dense medium, have been explicitly evaluated in previous works [17,18,16] without considering the weak 
magnetism corrections. We will discuss now the Hartree response functions. In this approach nucleon correlations are 
not taken into account. 

The corrections introduced by the inclusion of the weak magnetism term will appear in the polarization insertions 
and these can be separated into tensor, vector-tensor, and tensor-axial pieces which must be calculated in addition to 
the other well known vector, and axial-vector parts extensively quoted in previous works [17,18]. The polarizations 
induced by the weak magnetism in the Hartree approximation are: 

(12) 
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nS r = -i I tHj nCiip + q ) Y cm (13) 

= -i I ^Tr[Gi( P + q )(-ri 5 )G i {p)^^] (14) 
These are similar to the original vector- vector polarization [17], 

= -i J T ^( P + q ) r Gi( P ) r\ ■ (15) 

More information on these polarizations is in the appendix. At the Hartree level, the only modifications from the 
medium are Pauli blocking, and the density dependence of the nucleon effective mass. 

TABLE II. Coupling constants with g a = 1.260, sin 2 W = 0.231, fi p = 1.793 and fi n = - 1.913. 
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The cross section can be written as, 



cPcr _ G%E'ql 



d 2 n>dE> A^En B 



l-exp(-|J) (.17?, ,+R, ± /]/?,). (Hij 



with the plus sign for neutrinos, the minus sign for antineutrinos, and the kinematic factors, 



2EE' + \ql 



and, 



B = E + E'. (18) 

The three response functions are, 

* = E I (*? 2 + G ^ 2 ) I Im n -L + Im n« T ] + 2^ [Im H« L + Im H« T ] + if [Im H« + Im H« ] } , (19) 



i? 2 = E + °a ) Imn ™r + 2^ ImH« T + Ff IrnE^ - G^lmli^] , (20) 

Rs= 2 ( F i + F ijf) G A NlmnW. (21) 

i—p,n,e 

Hartree cross sections based on Eq. (16) with Eqs. (19,20,21) will be used to calculate neutrino mean free paths in 
Section III. 

B. RPA response 

We now discuss RPA corrections to the Hartree response functions from nucleon correlations. In a mean field 
approximation for the ground state of neutron rich matter, the meson fields are calculated self-consistently. This can 
be done with a Dyson equation that sums the tadpole diagrams. The linear response of this system can be obtained 
by attaching the weak interaction to all nucleon lines. This yields the RPA equations that sum the ring diagrams to 
calculate the nuclear response. Calculating cross sections in RPA yields neutrino opacities that are thermodynamically 
consistent with an equation of state calculated in a mean field approximation, see for example [16]. Therefore, the 
RPA opacities from this section are appropriate for supernova simulations involving the relativistic mean field TM1 
equation of state [13]. 

In the RPA, correction terms from correlations are added to the Hartree responses of Eqs. (19,20,21). 

R? PA = R 1 + AR L + AR T + AR aa (22) 

R PPA = i?2 + AR T + AR aa (23) 

Rf PA = R 3 + AR va (24) 



The cross section in RPA is, 

d 3 a G 2 F E'ql 



d 2 il'dE' 4n 3 En B 



l_ ex p(-|)] 1 (AR PPA + R PPA ±BR PPA ) . (25) 



The correction ARt from transversely polarized mesons and photons was discussed in ref. [17]. Here we add weak 
magnetism corrections, the p meson, and nonlinear meson self-couplings. 

AR T = Im|n^(l - D T n T )" 1 D T n TV1 /| (26) 
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The polarization IIt is a 3 x 3 matrix for the system of electrons, protons and neutrons, 



n 7 



W vvT 

o o n™, , 



(27) 



and the interaction matrix is, 



D 7 



-X 7 X 7 

Xi ~ Xv ~ Xp Xp ~~ Xv 

Xp ~ Xv ~Xv ~ Xp 



(28) 



The photon and meson propagators are, 



X 7 o ' 



(29) 



1.2 ' 
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where because of self-interactions the vector meson effective mass is, 



*2 2 i S> 2/ t/ \2 



(32) 



Here V is the vector meson mean field and the coupling constants are e 2 (e 2 /4-7r = a), g% 
The coupling of the system to weak currents is described by a three component vector U TW 



2 and g 2 p . 
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(33) 



and its transpose H-tw- This vector contains the weak magnetism corrections involving F\. 

The correction from axial currents Ai? aa is given by Eq. (50) of Ref. [17] with no corrections from weak magnetism. 

Ai? aa - q 2 Im{n* w (l - D T n T )" 1 D T n aW } (34) 

Here H a w is a three component vector describing the coupling to weak axial currents and n* w is its transpose, 



H a w — 



A^va 

G p A W va 

Gn TTTl 
A LL va 



(35) 



Note, we assume no additional corrections to Ai? aa from pion exchange or short range correlations that are often 
parameterized with a constant g' . These can significantly reduce the RPA cross sections at high density [11]. Therefore 
our results without g' represent a minimal calculation that is consistent with the TM1 equation of state. Clearly 
additional RPA interactions (such as g'), that do not contribute to the mean field equation of state, can be added 
without destroying the consistency. We discuss this further in section IV. 
The vector-axial correction AR va is, 



AR va = 2im{n* iy (i - D T n T )- 1 D T n TW } 



(36) 



Finally ARl describes corrections from longitudinally polarized mesons and photons. This correction is slightly 
more complicated because of mixing in the medium between scalar mesons and longitudinally polarized vector mesons. 
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This can be taken into account by increasing the dimension of the matrices from 3 x 3 to 4 x 4. The four rows describe 
vector coupling to electrons, protons and neutrons and a single scalar coupling. Because the weak interaction has 
no scalar coupling one can combine all of the scalar polarizations into a single entry that is the sum of proton and 
neutron contributions. This yields [17], 



ar l = im{n* LW (i - D L n L )- 1 D L n LW } 

Here the longitudinal polarization matrix is, 



(37) 
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The scalar polarization rig is the sum of proton and neutron contributions Us = U s + II™, with, 

d 4 p 



and the mixed scalar- vector polarization is, 



(M 



: Tr[Gi(p + q)Gi(p)], 



n 



(39) 



(40) 



The interaction matrix is, 



D 7 



with the scalar meson propagator, 



The scalar effective mass is, 
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(43) 



Finally, the longitudinal coupling of the system to the weak current is a four component vector, 



n 



LW 



FfKvL 

FfKvL + F 2 p K tL 



(44) 



and its transpose is n^^. Note, the second line uses results for the mixed scalar-tensor polarization, see the appendix. 
It is a simple matter to numerically evaluate the corrections ARl, ARt, AR aa and AR va to calculate the RPA cross 
section. 



C. v Opacities 

We now calculate the transport mean free path for /x or r neutrinos or antineutrinos in dense neutron rich matter. 
We assume the opacity is dominated by neutral current scattering from protons, neutrons and electrons. Well inside 
the neutrino sphere, the mean free path becomes much smaller than the size of the system, and we expect a simple 
diffusion approximation to be valid. Using the diffusion approximation, the energy flux F can be written 
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(45) 



where U is the energy density of the neutrinos. 

The transport cross section for a neutrino of energy E scattering through an angle 9 is, 



a ' {E) = I dE> / d "'^S^ (1_COsW) ' (46) 



and the mean free path is simply, 



A < £ > - (47) 



Note, when evaluating the integration over the outgoing neutrino energy E' , the energy transfer E — E' takes both 
positive and negative values. 

The mean free path averaged over the incoming neutrino energy spectrum is, 

£ X(E)E\f(E)dE 
K A >_ J™E*f{l3)dE ' (48) 

For simplicity we assume f(E) w cT E l T at temperature T. The mean free path for a neutrino of average energy, 
A(< E >) with < E >= ttT, and the mean free path averaged over energy < A > are in general somewhat different. 
To quantify the differences between antineutrinos and neutrinos, we define a ratio, 

R = <I> ; <A> (49) 
< A > 

with < A > the antineutrino and < A > the neutrino mean free path. This difference comes mostly from weak 
magnetism corrections. 



III. RESULTS 



We now present results for mu and tau neutrinos and antineutrinos scattering in neutron rich matter. At low 
densities and temperatures there is a transition from a uniform liquid to a nonuniform solid phase. This solid phase 
is expected to make up the crust of a neutron star. The instability of the uniform system to density oscillations 
shows up as an enhancement in the RPA cross section for small values of the energy transfer qo — E — E' where E 
is the initial and E 1 the final neutrino energy. Note, all of our calculations assume a uniform system. Therefore they 
are not directly valid in the nonuniform phase. However, the enhancement of the cross section we find in an RPA 
approximation could be a first hint of the perhaps greatly enhanced neutrino opacities we expect for the nonuniform 
phase, see below. 

The differential cross section da/dfldE' for a 20 MeV neutrino scattering at 90 degrees from neutron rich matter 
is shown in Figures (1,2,3) at three low densities. The electron fraction is Y e = 0.3 and the temperature is T = 5 
MeV. At a density of p = 0.01 fm~ 3 the uniform phase is unstable. This shows up as a large peak in the RPA cross 
section near q = 0. This peak decreases as the density increases until it almost disappears by p = 0.1 fm~ 3 where 
the uniform phase is stable. 

Explicit calculations of opacities for the nonuniform phase would be very useful and are an important area for future 
work. We expect the opacity to be greatly enhanced because (1) neutrinos can scatter coherently off of the density 
changes. For example Reddy et al find greatly increased opacities for a nonuniform Kaon condensed phase [23]. We 
note that the neutrons in neutron rich matter carry large weak charges and the neutrino wavelength is comparable 
to the size of the density fluctuations. (2) There should be many collective modes at low excitation energies because 
there may be many different "pasta" configurations and shapes that all have very similar energies. The first hint of 
these modes shows up in RPA calculations as the excitation energy of collective density oscillations drops to zero. 
Neutrinos may be able to efficiently excite these modes. Furthermore, neutral current excitation of these modes may 
allow mu and tau neutrinos to more efficiently exchange energy with matter and remain in thermal equilibrium to 
lower densities and temperatures. 

Figure (4) shows the transport mean free path in a Hartree approximation for 9^ and v T versus density. The energy 
averaged mean free path Eq. (48) is seen to be somewhat smaller than the mean free path evaluated at an average 
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energy E = irT. This could be because Pauli blocking changes the energy dependence of X(E) from 1/E 2 . Note, in 
all of the following figures we plot energy averaged quantities. 

Figure (5) shows the transport mean free path for scattering from neutron rich matter in beta equilibrium at a 
temperature of 5 MeV. Note that the mean free path actually increases with density between ns ~ 0.2 and « 0.5 
fm~ 3 because of the rapid decrease of M* with density. Weak magnetism causes the mean free path for neutrinos to 
be systematically smaller than that for antineutrinos. At high densities, RPA correlations increase the mean free path 
somewhat. Our calculations use the simplest RPA interaction that is consistent with the TM1 equation of state. We 
do not include residual interactions from tt exchange or the effects of short range correlations. These may increase 
the RPA mean free path further at high densities [11]. Unfortunately, the parameter g' which was used in [11] to 
include short range correlations may have some density dependence. This should be investigated in future work. 
Nevertheless, we expect the significant differences between neutrino and antineutrino mean free paths to persist in 
calculations using more sophisticated interactions. At lower densities, the RPA mean free path becomes shorter than 
that in a Hartree approximation because of coupling to density fluctuations related to the transition to a nonuniform 
phase. This was previously discussed. 

Increasing the temperature to 10 MeV, as shown in Figure (6), increases the average neutrino energy. Weak 
magnetism effects grow with energy so the difference between neutrino and antineutrino mean free paths is now larger 
than at T=5 MeV. This difference continues to grow as the temperature is increased further to 25 MeV in Figure (7) 
and 50 MeV in Figure (8). At these higher temperatures there is no longer a significant reduction in the RPA relative 
to Hartree mean free paths at low densities. This is because the system is now far from any phase transition. 

Figure (9) shows the temperature dependence of the mean free paths and Figure (10) shows the fractional difference 
R between antineutrino and neutrino mean free paths versus temperature, see Eq. (49). This difference is large and 
grows rapidly with temperature. At high temperatures the antineutrinos have much larger mean free paths. Our 
RPA calculation has a weak density dependence. This suggests that the simple weak magnetism correction factors 
calculated in Ref. [4] at zero density may be a reasonable first approximation to results including correlations in dense 
matter. The large ratio R may lead to and v T spectra that are somewhat colder than spectra for and 9 T . 

IV. CONCLUSIONS 

Advances in neutrino transport allow realistic neutrino interactions to be incorporated into supernova simulations. 
This is important because of the very large energy radiated in neutrinos. We believe a minimum requirement for 
realistic neutrino opacities is that: (a) they agree, at low densities, with free space neutrino-nucleon cross sections 
that include weak magnetism and recoil corrections, (b) Neutrino opacities should be consistent with the equation of 
state. To our knowledge, no present simulations have satisfied these basic requirements. 

In this paper we have added weak magnetism corrections to relativistic RPA calculations. Our results agree at 
low density with accurate neutrino-nucleon cross sections including the fundamental differences between neutrino and 
antineutrino interactions. Furthermore, our opacities are consistent with the relativistic mean field equation of state 
based on the TM1 interaction. 

We find antineutrino mean free paths that are considerably larger than those for neutrinos. This is especially true 
at high temperatures where neutrino energies are large. This difference, between v and v mean free paths, depends 
only weakly on density and could produce and v T spectra that are somewhat colder than those of and v T . 

We find RPA correlations somewhat increase mean free paths at high densities. However, we do not include any 
residual interaction from pion exchange or short range correlations. These; could increase the mean free path further. 
Therefore our results represent a minimum calculation that is consistent with the equation of state. Future work 
studying the role of short range correlations and the density dependence of the effective interaction would be very 
useful. 

At low density, correlations reduce the mean free path because of excitation of low energy collective density os- 
cillations. These are related to the eventual transition of neutron rich matter to a nonuniform phase. Our RPA 
calculations can only describe very small amplitude density fluctuations. Nevertheless we expect neutrinos will couple 
very effectively to the large amplitude density fluctuations present in the nonuniform phase. This could lead to a 
large decrease in the neutrino mean free path and change the emitted neutrino spectrum. We believe, it is important 
to explicitly calculate neutrino opacities for the nonuniform phase. 
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V. APPENDIX: POLARIZATION INSERTIONS WITH TENSOR COUPLINGS 

In this section we provide relations for polarization insertions with tensor couplings. These compliment previously 
published results for polarizations with vector couplings. See for example [18] at finite temperature and [17,24] at 
zero temperature. We use the axial-axial polarization, 

K: = ~if (04 Tr t G (P + Qh^GiPh^], (50) 

that can be written, 

n£ = n{£ + <rn A) (5i) 

which defines 11^. Using this we have, 

U vtT = -U vtL = ^-Il A . (52) 
For tensor-tensor correlations we will only need the imaginary parts. 

% 



2 .. „2 



ImHttT = ^-rrlm 

4M 2 



ImI1 ^ = -i^ Im 



IW - (1 + ^ )Ua (53) 



„2 



Tl vvL + (l + ^)Tl A \ (54) 



The vector-axial polarization, 



can be written, 



Kl = -^ a Vn™ . (56) 

The tensor- axial polarization is simply related to ILj , 

M* 

Tlta = w Tl va . (57) 

This has been used in Eq. (21). Finally, the mixed scalar-tensor polarization is simply related to the scalar- vector 
polarization, 

n ^ = mw n -- (58) 

This has been incorporated into Eq. (44). 
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FIG. 1. Differential cross section dcr/d£ldE' versus excitation energy qo = E v — E' v for a baryon density p — 0.01/m~ 3 , a 
neutrino energy E v — 20 MeV, a scattering angle of 90 degrees, an electron fraction Y e = 0.3 and a temperature of 5MeV. The 
solid curve is in a Hartree approximation while the dashed curve includes RPA correlations. 
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FIG. 3. As per figure 1 at a baryon density of p = 0.1 fm 




FIG. 4. Transport mean free path in meters for or v T at a temperature of T=10 MeV as a function of baryon density for 
neutron rich matter in beta equilibrium. The solid curve is for a neutrino energy of E v = 7vT — 31.4 MeV while the dotted 
curve is the mean free path averaged over neutrino energy. Both calculations are in a hartree approximation. 
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FIG. 5. Transport mean free path in meters averaged over neutrino energy for neutron rich matter in beta equilibrium at 
a temperature of T = 5 MeV. The solid curves are for and v T while the dotted curves are for or v T . At high density, 
the upper curves are RPA and the lower curves Hartree calculations. Note, that the RPA mean free paths are smaller at low 
density than the Hartree ones. 




FIG. 6. Transport mean free path in meters averaged over neutrino energy for neutron rich matter in beta equilibrium at 
a temperature of T = 10 MeV. The solid curves are for and v T while the dotted curves are for or v T - At high density, 
the upper curves are RPA and the lower curves Hartree calculations. Note, that the RPA mean free paths are smaller at low 
density than the Hartree ones. 
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FIG. 7. Transport mean free path in meters averaged over neutrino energy for neutron rich matter in beta equilibrium at a 
temperature of T — 25 MeV. The solid curves are for and v T while the dotted curves are for or v T . At high density, the 
upper curves are RPA and the lower curves Hartree calculations. 
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FIG. 8. Transport mean free path in meters averaged over neutrino energy for neutron rich matter in beta equilibrium at a 
temperature of T = 50 MeV. The solid curves are for and v T while the dotted curves are for or v T . At high density, the 
upper curves are RPA and the lower curves Hartree calculations. 
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FIG. 9. Transport mean free path in meters averaged over neutrino energy for neutron rich matter in beta equilibrium versus 
temperature in an RPA approximation. The baryon density is 0.05 /m~ 3 for the solid curves, 0.15 fm~ 3 for the dashed and 
0.53 fm~ 3 for the dotted curves. The upper curves are for i/ M and v T while the lower curves are for and v T . 




FIG. 10. Fractional difference between antineutrino and neutrino mean free paths averaged over energy versus temperature 
in an RPA approximation. See Eq. (49). The baryon density is 0.05 fm~ 3 for the solid, 0.15 fm~ 3 for the dashed and 0.53 
fm~ 3 for the dotted curve. 
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